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A two-stroke, externally heated valve engine (EHVE) with a heater, a cooler and two blowers is simulated. 
The engine is entirely different from a typical Stirling engine. The pressure ratio p max /p min of its cycle is 
higher, but the engine volume and the mean value of the heat exchanger wall temperatures are the same. 
The power and efficiency of the EHVE and Stirling engines under the same maximum pressures are 
compared. The results show that the EHVE engine reaches almost the same level of performance as 
the Stirling engine, while using only available atmospheric air, rather than helium. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

In this paper, a new type of externally heated valve engine 
(EHVE), fully described in [1 ], is examined and quantitative results 
of its performance are presented. Schematic views of the engine 
are shown in Figs. 1 and 2. 

Stirling engines are well known externally heated devices and 
remain a subject of current research [2-6], Such engines have sim¬ 
ilarities to the EHVE, including: external heating; high thermal effi¬ 
ciency; multi-fuel capability; low noise generation; and low 
pollution when a combustion process is used as the heat source. 
A general feature of the EHVE and Stirling engines is the ability 
to use a range of heat sources. When a combustion process is used 
for this purpose, the burning arrangement may be optimised for a 
particular fuel. However, heat sources are not limited to combus¬ 
tion processes only, but can involve any type of energy capable 
of reaching the required temperature. 

There are many papers showing experimental results of differ¬ 
ent Stirling (or similar externally heated) engines designed to use 
various alternative sources of energy, see [7-24], although not all 
were intended to reach high performance levels. Applications 
range from demonstration devices such as [23], through general 
purpose machines [7,21], solar heated [13], marine use [17,26], 
to potential outer space energy conversion [25], The reader can 
find a comprehensive review of such solutions in [26-29] and a 
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summary of the performance of a number of existing machines 
in [4[. 

Numerical simulations of the EHVE operation are compared 
with the experimental results presented for the Stirling engine de¬ 
scribed in [2,3] (1987-1988). The engine has a compact mechanical 
form and comparatively high power and efficiency at a total 
volume of 588 cm 3 enclosed in four double-action cylinders with 
four piston rods. The working gas is helium. Later Stirling engines 
described in [4-6] (2001-2003) have analogous or even lower 
effectiveness at similar design parameters. An initial comparison 
of the EHVE and an example of the Stirling engine has been made 
already in [30], but taking into account the later results developed 
in [31 ], such evaluation could be too optimistic for the EHVE. That 
evaluation is expanded and completed here. 


2. Design aspects of the EHVE 

The EHVE in the earlier form as shown in [30] exists as a pro¬ 
totype version ([31 , Fig. 3). The experimental results have al¬ 
ready proven the concept. In the prototype, an electric heating 
device was applied. However, independently of its type, it re¬ 
sulted in insufficient heat being delivered during the period when 
the engine valves were closed. An intensive flow of air is required 
in the heater throughout the complete cycle. Therefore, two cen¬ 
trifugal blowers Eh and B 2 were added to the model of the EHVE 
in order to improve the heat exchange rate. Their power require¬ 
ments are minimal but raise the heat exchange level when the 
engine valves are closed. The version of the EHVE discussed here 
is equipped with a single heater and a single cooler. The volumes 
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Nomenclature 



A 

surface area (m 2 ) 

>1 

efficiency 

D, d 

diameters (m) 

K 

isentropic expansion index 

hsBi 

enthalpy isentropic increase of blowers (J/kg) 

CD 

angular velocity of the crankshaft (rad/s) 

L 

mechanical work (J) 

Subscripts B u B 2 blowers 

1 

length of tubes (m) 

C 

compressor 

rhj 

mass flow rate (kg/s) 

Cl 

cooler 

P 

pressure (N/m 2 , bar) 

Ch, cu 

collectors of the cooler 

P 

power (W, kW) 

E 

expander 

n 

number of tubes 

f 

friction 

Q. 

heat (J) 

H 

heater 

Qa 

heat flow in the surface area (W) 

h,,h 3 

collector of the heater 

R 

specific gas constant (J/kg I<) 

i 

number of the blowers (=1, 2) 

Re 

Reynolds number 

j 

number of the valve (=1, 2, 3, 4) 

S 

entropy (J/kg I<) 

k 

opening (k = 0), closing (k = c) 

T 

temperature (I<) 

max 

maximum 

V 

volume (m 3 ) 

PH 

pre-heater 

a 

crankshaft angle (rad, °) 

w 

wall 

® j,k 

angle of valves (rad) 





Fig. 1. Scheme of the transversal cross-section of the engine E - expander, C - 
compressor, H - heater. Cl - cooler, Bi, B 2 - blowers, 1, 2, 3, 4 - valves, Hi, H 3 - 
heater collectors, Cl 2 , Cl 4 - cooler collectors, PH - pre-heater. 

of each heat exchanger in the EHVE are at least 3-5 times larger 
than the volumes of the cylinders. The inlet and outlet collectors 
of the heater and the cooler are shown in Fig. 1. Their volumes 
may play a role of settling chambers for the working gas flow 
(equalizing pressure and temperature values in the exchangers 
in time). 

An important advantage of this configuration is use of the most 
easily available working gas-atmospheric air. Also, it is possible to 
utilise any source of heat able to deliver an appropriate tempera¬ 
ture to the heater, that is, from solar energy, through combustion 
of other media, to nuclear power. This makes the EHVE a suitable 
solution when available resources of oil and gas run out. The high 
temperature of the heating gas leaving the EHVE allows for the 
installation of an additional heat exchanger, a pre-heater (PH), 
improving thus efficiency, see Fig. 1 and Section 3. The earlier 
EHVE scheme, shown in [1], did not include such a device. 




Fig. 3. General view of the EHVE prototype. Engine block in the middle, heat 
exchangers at sides - two electric heaters have been used in the experiments 
performed in 2001, see [4]. 
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The EHVE used in the comparison is a two-stroke type device. It 
is equipped with two cylinders, each of which plays a different role 
(Figs. 1 and 2). The expander E and the compressor C are situated 
between two external heat exchangers, the heater El and the cooler 
Cl, with the phase angle between them equal to 180°. Both of the 
heat exchangers, H and Cl, are additionally equipped with their 
own blowers Bt and B 2 raising the rate of heat exchange. The gas 
flow inside the exchangers is assumed to be turbulent even when 
the engine valves are closed. The engine works in a closed cycle. Its 
valves are operated by means of a typical camshaft. Additionally, 
the EEIVE features a conventional crankshaft and a conventional 
oil lubrication system, which is an advantage for manufacturing. 
As in the internal combustion (IC) and existing Stirling devices, 
one can expect an increase of performance in a multi-cylinder 
arrangement of the engine [9], 

The pressure ratio p max /Pmm °f the EHVE and the Stirling engines 
considered is different. For the EHVE, it is greater than 5, while for 
both of the Stirling engines considered, it is less than 2. For this 
reason, comparisons cannot be performed with equivalent mini¬ 
mum pressures. For the Stirling engines, the minimum pressures 
are equal to 60, 82 and 110 bar, but in the EHVE such minimum 
pressures would imply unrealistically high values of the respective 
maximum pressures, which may be not acceptable from an engi¬ 
neering point of view. The new form of the EHVE with blowers is 
therefore compared with the Stirling engines using the same max¬ 
imum working pressure. 

3. Simulations of the EHVE 

The model of the EHVE and its mathematical realisation is de¬ 
scribed in detail in [1], We repeat only the important elements 
here. These are time-dependent relations for governing the engine 
action, in which dependence on spatial variables has been omitted, 
so the model is described and simulated in the time domain only. 
In both heat exchangers, the spatial changes are important be¬ 
cause, e.g., the temperature of the working air inside the exchang¬ 
ers in fact varies on the way from the inlet to outlet collectors. A 
simpler approach to the problem is also possible and it has been 
applied in [1], 

The heat exchangers can be described as follows. The heater and 
the cooler are counter-flow heat exchangers, where the tempera¬ 
ture gradient remains nearly constant along their length. The heat 
exchangers are treated as time-dependent only, assuming that all 
important changes are the same along the overall length of the 
exchangers and considering each of them in a single plane only. 
The plane is situated at the outlet of the exchanger tubes. The ap¬ 
proach used here assumes the application of maximum wall tem¬ 
peratures of the heater tubes and minimum temperatures of the 
cooler wall tubes. This method of analysis of heat exchangers is 
sufficient to demonstrate the advantages resulting from use of 
such blowers. Such a treatment of the parameters is an approxima¬ 
tion, but it allows the use of a moderate number of equations in the 
time-dependent model. 

The blower in the loop of hot and cool air circulation ensures 
that the Reynolds number calculated for the exchanger tubes is 
not less than 10 4 , when both valves 1 and 3 remain closed. 

The work of the expander and the compressor is calculated for 
an angle 0° ^ a ^ 360°, in the same way as for the 2-stroke engine: 

,*=360° dV E , 

Le= .L P£(a) w da ’ (1) 

,*=360° 

u= L pc(a) w da - (2) 

The pressures in the above formulas are determined according 
to [1], taking into account losses in the valves. 


The work of the engine is the sum of: 

L = L e + Lc — Lbi — £b2 — Lf — L sup , 
where the values of work of the blowers are (i = 1,2): 
2mh Bi h S Bi 

L Bi ~ -■ 


(3) 

(4) 


L f is the work due to friction in the machine, and, L sup is the work 
lost through all additional devices attached to the EHVE. Frictional 
losses in modern IC engines are observed to be between 3% and 
8% of total horsepower. As the EHVE can use almost the same oil 
lubrication system as an IC engine, similar frictional loss values 
were applied. Its state-of-the-art design can reduce amount of the 
losses considerably. 

The symbols used are explained in the nomenclature, and an 
enthalpy isentropic increase of blowers is calculated from: 

hsB,=^RT bB ,(nf-L). (5) 

A complete power output of the externally heated 2-stroke en¬ 
gine is determined as: 


P =S L 

The heat delivered inside the heater (following the theoretical 
model described in [1]) is: 


Qh = 


/»a=3 

J a=0° 


Qah 


da 
a>' 


(7) 


The efficiencies of the engine may be defined as: 




( 8 ) 


The calculation of the value of Q H is made without an additional 
heat exchanger. However, the temperature of the gas heating the 
heater tubes at its outlet (Fig. 1) would be very high, so a pre-hea¬ 
ter may be applied as described in [2,3], The pre-heater being an 
additional heat exchanger raises the inlet temperature of the heat¬ 
ing gas and causes considerable savings in heat requirements 
(decreasing the denominator value in Eq. (8) from Q H to Q H _ PH ). 
The efficiency of the EHVE should increase significantly in this case. 


4. Comparison of the air EHVE to helium Stirling engines 

The numerical simulation of the EHVE then requires a solution 
of 19 non-linear ODE’s with appropriate switching conditions, sup¬ 
plemented by a set of nonlinear algebraic equations describing the 
behaviour of elements of the model. Stationary, periodic solutions 
to the engine operation at a given, constant angular velocity of the 
engine are presented. The simulation program uses a C++ language 
library of procedures and solvers (Runge-Kutta fourth-order meth¬ 
od) [32]. All simulations started from identical, fixed initial condi¬ 
tions and continued until the solution reached the point where the 
integrated values differed by no more than 0.01% from those ob¬ 
tained at the end of the previous cycle. A good agreement of the 
integrated values with the experimental results for the prototype 
engine was found, as published earlier in [31]. Typically, a stable 
solution was reached after 150 cycles of simulated crankshaft 
revolutions. 

The calculations were performed for the ratio of volumes 
V C /V E = 0.75, which was found to be the optimal value for the case 
when V H ss 5 V E , see [1 ]. The EHVE has the same expander volume 
V E as the Stirling engine [2,3], which is a double-action system. 
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Other essential data, different to those presented in [1], are as 
follows: 

d E = 0.0900 m, s E = 0.0925 m, h 0E = 0.0025 m, 

d c = 0.0900 m, s c = 0.0694 m, h oc = 0.0017 m 

- representing dimensions in the volumes of E and C. 

Dimensions of the heat exchangers (d - pipe diameter, / - pipe 
length, n - number of pipes, and, V - the exchanger volume) have 
been assumed as: 

d H = 0.01 m, l H = 2.00 m, n H = 15, V m = V H3 = 0.001 m 3 , 

d a = 0.01 m, la — 2.00 m. n a = 15, V C a = V C i 4 = 0.001 m 3 . 

The value of A p H3 depends on the design of the collectors and 
the heat exchangers. Here, all values of the Ap’s have been fixed 
at 10 kPa. The mass flow rates of the blowers have been estimated 
at ni B i = 0.1 and 0.2 kg/s. Other important blower parameters are: 
input to output pressure ratio /7 B , = 1.01, and efficiency t] Bi = 0.7. 

The work of friction during the cycle is estimated as (derived 
and adjusted to the EHVE from [33-35]): 

L f = 4 0+08 £- 

The power for all auxiliary devices - L sup , is supplied by an electric 
motor rotating at 3000 rpm, independently of the EHVE rotation, 
and is also taken into account. For this purpose, the calculation as¬ 
sumes a power consumption of 1000W. £4 

The temperature of the EHVE heater walls is assumed to in¬ 
crease with the engine rotational speed from 1100 to 1250 K and 
it is related to the amount of hot heating gas delivered, see Table 1. 

This temperature reflects the limited amount of data given in [2,3] 
for the Stirling engines. The minimum temperature of the cooler 
wall is constant and equal to T wa = 298 K based on standard 
parameters of the cooling water supply and data from [2,3], The 
coefficients of heat conduction, dynamic viscosity and specific 
heats are taken from tables of physical properties of air. 

The simulation program, based on the model of the EHVE de¬ 
scribed in [1] has been used for calculations for the parameters 
reflecting the data available for the two Stirling engines presented 
in [2,3] at selected rotational speeds from 500 to 2000 rpm. 

In Fig. 4, a difference in the area enclosed by the p-V curves can 
be interpreted as the useful work provided by the EHVE. Fig. 5 
illustrates a T-S diagram of the EHVE and a short description of 
the operation cycles. The cycle shown in Fig. 5 is a new type of 
the heat exchange cycle. The statement that it is similar to the 
Joule cycle, which was given in [1], should be treated as a 
simplification. 

Use of a pre-heater does of course reduce the energy require¬ 
ments of the heater. All details presented in [2], (earlier than 
[3]), and even the pre-heater described there are different when 
compared to [2], Therefore, it has been assumed Q H -ph = 0.80Q H 
for [3] and Q H -ph = 0.85Q H for [2], These assumptions are rather 
modest and can be easily met in real installations. The efficiency 
of the EHVE in this case is obviously higher than in the case with¬ 
out a pre-heater and is presented in Figs. 6-8. It is important that 
the Stirling engines against which the EHVE is compared make use 
of pre-heaters [2,3], The heat delivered to the EHVE reflects the 

Table 1 

Heater wall maximum temperatures of the EHVE as a function of rpm. The heater wall 
mean temperature is the same for both the air EHVE and helium Stirling engines [2,3] 
and equal to 1023 K. 

n E HV e (rpm) 500 1000 1500 2000 

T wH max (I<) 1100 1150 1200 1250 



F[m 3 ] 


Fig. 4. p—V diagram for the EHVE at 1500 rpm, p max = 140 bar and = 0.1 kg/s. 
Data in Section 3. 



Fig. 5. Approximate T—S diagram for the EHVE at 1500 rpm, p max = 140 bar and 
m Bl = 0.1 kg/s. 1-2: isentropic compression inside the closed C volume, 2-3: 
compression by C, heating by H and mixing in the collector H 3 by the blower B lt 3- 
4: heating inside the H volume, 4-5: expansion in E with heating at opened valve 1 
until its closing, 5-6: isentropic expansion inside the closed E volume, 6-7: flow to 
the collector Cl 2 , expansion and mixing with the collector Cl 2 by the blower B 2 ,7-1: 
cooling inside the Cl volume. 



Fig. 6. Comparison of the EHVE and AIS1N Stirling engines at the maximum working 
pressure p max = 102 bar. Power values presented with empty markers and effi¬ 
ciency with the filled ones. Circles on solid lines represent the EHVE with its 
blowers working at m Bi = 0.1 kg/s, top-triangles in grey, dashed line with blowers 
at m B i = 0.2 kg/s, whereas star symbols stand for the AISIN engine [3]. 
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rpm 

Fig. 7. Comparison of the EHVE and AISIN Stirling engines at the maximum working 
pressure p max = 140bar. Power values presented with empty markers and effi¬ 
ciency with the filled ones. Circles on solid lines represent the EHVE with its 
blowers working at m Bi = 0.1 kg/s, top-triangles in grey, dashed line with blowers 
at m Bi = 0.2 kg/s, whereas star symbols stand for the AISIN engine [2]. 



Fig. 8. Comparison of the EHVE and AISIN Stirling engines at the maximum working 
pressure p max = 190bar. Power values presented with empty markers and effi¬ 
ciency with the filled ones. Circles on solid lines represent the EHVE with its 
blowers working at m Bi = 0.1 kg/s, top-triangles in grey, dashed line with blowers 
at m Bi = 0.2 kg/s, whereas star symbols stand for the AISIN engine [2]. 


unimproved pre-heater and is then assumed to be equal to 
Qh-ph = 0.85Q h . 

The comparisons presented below have been made between the 
numerical results obtained here and the experimental results re¬ 
ferred to the helium Stirling engine described in [2,3]. The results 
are presented in Figs. 6-8 for the same engine volume, rotational 
frequency and maximum pressures. The maximum pressure values 
of 102,140 and 190 bar for the Stirling engine cycle have been de¬ 
rived from the data given in [2,3]. Taking into account the results 
presented in [3], the heat exchangers, including a pre-heater are 
improved with efficiency values about 4% higher than those quoted 
in [2], 

The blowers B 3 and B 2 working at a flow rate of m B , = 0.1 kg/s 
are necessary for efficient operation of the EHVE. A further increase 
in the blower flow rate is not practical, see (Figs. 6-8). The reason is 
an increase in the power consumed by both blowers. The compar¬ 
ison also shows that the EHVE for all values of p max at 
m B i = 0.1 kg/s can produce slightly more power than the equiva¬ 
lent helium Stirling engine. It must be noted that this is achieved 
with air as a working gas rather than helium. 


A comparison of the efficiency levels of both engines is more 
complicated. The maximum efficiency of the Stirling type shows 
a distinct peak at 1000 rpm, whereas the efficiency plot for the 
EHVE is rather flatter. In Fig. 6, where heat exchangers (basic type 
in the case of Stirling engines) are almost the same, it is seen the 
efficiency of the EHVE is nearly equal. In Figs. 7 and 8, where heat 
exchangers in the Stirling engine are improved, the EHVE efficiency 
is lower on average by about 4%. 

The presence of valves in the EHVE also causes some disadvan¬ 
tages, which are as follows: 

(1) a short period of opening of valves 1 and 3, as in [1], 
Aoq = 90°, and Aa 3 = 80°. 

(2) the blower B 3 works at high temperature and pressure con¬ 
ditions (Fig. 1). 

The period of the valve opening phase is a challenging but fea¬ 
sible engineering task. The necessary work related to valve timing 
is still in progress [36]. The heater, not discussed here in detail, is 
also a demanding component of the EHVE. High working pressure 
and temperature are the critical factors in its design. Continuous 
improvements of heaters used in the Stirling engines however 
allow us to anticipate their availability for the EHVE without any 
major changes. Some examples of heaters and heat exchange prob¬ 
lems are presented in [37-40]. 


5. Conclusions 

A two-stroke externally heated valve engine called the EHVE, 
entirely different from the well known Stirling type, has been 
investigated by numerical simulations. Two counter-flow heat 
exchangers, a heater and a cooler, and a pre-heater are employed 
in the configuration. These have been supplemented by two blow¬ 
ers Bj and B 2 in order to increase the rate of heat exchange. 

The results have been compared with the experimental data ob¬ 
tained for very good examples of helium Stirling engines [2,3] for 
similar maximum pressures of their respective cycles. The same 
engine capacities, heater wall temperatures and working frequen¬ 
cies have been assumed, see Figs. 6-8. The differences in pressure 
ratios p max /p min of the compared engines being about 5:2 are 
mainly caused by the presence of valves and design of dead space. 
The Stirling engine has no valves and large dead space, contrary to 
the configuration of the EHVE. 

The EHVE presents an improvement in the area of externally 
heated engines, especially over the Stirling designs in terms of 
power, losing only some efficiency under similar working condi¬ 
tions. The most important advantage of the new design is such that 
these results have been achieved using atmospheric air rather than 
helium. We estimate that using helium, the EHVE would improve 
its power and efficiency by further 30% compared to the results 
presented. The EHVE design also has the advantage of less compli¬ 
cated technology at the manufacturing stage than that required for 
Stirling engines. The use of modern manufacturing solutions, as 
applied to typical internal combustion engines, would allow a con¬ 
struction of the EHVE. 

One can speculate that EHVEs might be used on small nuclear 
powered submarines as they are as quiet as Stirling engines. Strong 
sunlight may deliver sufficient energy to run the EHVE. The power 
provided by this engine could be further increased through the use 
of a tracking device to maximise the solar energy harvested. Con¬ 
sidering the above mentioned advantages, it can be expected that 
this new engine will raise interest among the world energy com¬ 
munity. In other words, the proposed engine combines the best 
of two worlds - the well-known technology of manufacturing of 
internal combustion engines and an environment friendly external 
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burning process (among other available sources of heat) - and can 
stimulate a wide range of applications of externally heated engines 
in the near future at a reasonable cost. 

A lack of any valves is an important advantage of the Stirling en¬ 
gine, but we need to emphasise the advantages of the EHVE: 

(1) the working gas is air, whereas the Stirling engines consid¬ 
ered have to use helium to get a comparable level of power 
and efficiency; 

(2) the EHVE uses a less complicated, conventional crankshaft 
system, see Fig. 2, met in most modern internal combustion 
engines; 

(3) lack of a heat regenerator (additional heat exchanger), one of 
the essential parts of Stirling engines, in the EHVE makes the 
heat exchange system simpler; 

(4) another advantage of not using a regenerator (usually made 
of steel wires) is such that the EHVE can be lubricated with 
the same oil system as in internal combustion engines. 
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